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Abstract 

We study the classical dynamics of the collinear positron-hydrogen scattering system below the 
three-body breakup threshold. Observing the chaotic behavior of scattering time signals, we in¬ 
troduce a code system appropriate to a coarse grained description of the dynamics. And, for the 
purpose of systematic analysis of the phase space structure, a surface of section is introduced being 
chosen to match the code system. Partition of the surface of section leads us to a surprising conjec¬ 
ture that the topological structure of the phase space of the system is invariant under exchange of 
the dynamical variables of proton with those of positron. It is also found that there is a finite set 
of forbidden patterns of symbol sequences. And the shortest periodic orbit is found to be stable, 
around which invariant tori form an island of stability in the chaotic sea. Finally we discuss a 
possible quantum manifestation of the classical phase space structure relevant to resonances in 
scattering cross sections. 
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I. INTRODUCTION 


The three-body problem of Coulombic systems is one of the most fundamental problem 
of theoretical physics. A large number of papers has been published which deals with the 
classical mechanics as well as the quantum mechanics of electron-hydrogenic atom systems, 
that is two-electron atoms dEl. The classical motion of three bodies in a two-electron atom 
above the three-body breakup threshold is known to be regular although it is nonintegrable 
[1], and the behavior of electron-impact ionization of hydrogenic atoms near the breakup 
threshold has been treated long time [HE]. In contrast to the dynamics above the threshold, 
the classical motion below the threshold is no longer regular; the classical phase space consists 
of not only reguar but also chaotic parts [7j. The chaotic dynamics below the threshold 
produces inhnitely many resonances in photoelectron spectra of two-electron atoms, and 
the fluctuating parts of the spectra owing to the resonances show a series of seemingly 
meaningless irregular signals [8]. However, by the help of modern techniques of semiclassical 
mechanics, the fluctuating parts can be characterized by the actions and stabilities of the 
inhnite number of closed triple collision orbits (CTCO) [31 E]. 

Positron-hydrogen scattering system is another three-body Coulomb problem, which has 
recently attracted a growing interest as powerful positron beams have developed since the 
discovery of low energy positron beams [TOHTS] . The behavior positron-impact ionization 
of hydrogenic atoms above the breakup threshold, whose corresponding classical dynamics 
is regular as in the electron impact ionization, has been examined by both classical and 
quantum mechanical methods [T3Hn]. Considering there are a large amount of papers for 
the classical dynamics of two-electron atoms for wide range of energies including below 
three-body breakup threshold, however, it is noteworthy that only a limited number of 
papers deal with classical dynamics of positron scattering by a hydrogen target below the 
threshold, where the classical motion might be no longer regular and thus fluctuations in 
quantum scattering cross sections might occur as in the case of two-electron atoms. 

Stimulated by a prominent role of classical calculations in success of semiclassical inter¬ 
pretations and predictions on quantum mechanical properties of two-electron atoms [7|, a 
study was conducted of classical dynamics of positron-hydrogen scattering for wide range 
of energies including below as well as above the threshold [18]. However, in that work, the 
classical mechanics was merely used as a calculational method for Monte Carlo simulations 
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rather than being analyzed for the pnrpose of characterization of the resonance strnctnres 
as was done in two-electron atoms. A systematic analysis of the classical phase space of the 
positron-hydrogen scattering system has not been given yet for the energies relevant to the 
resonances, that is for energies below the breaknp threshold 

Resonances lying close to the threshold are assnmed to be associated with CTCOs as 
demonstrated in papers deals with photoelectron spectra of two-electron atoms [alElEQ]. 
By virtne of the fact that CTCOs are completely inclnded in the collinear invariant snbspace, 
stndy of classical dynamics in the snbspace would give a great insight into the resonances 
in positron-hydrogen scattering in the real three-dimensional world. 

In this paper we will report on our investigation of the classical dynamics of positron 
scattering by a hydrogen target in the collinear conhguration for energies below the 

three-body breakup threshold. The paper is organized as follows. In Sec. II the Hamiltonian 
is introduced, and a transformation is also introduced to regularize Coulomb singularity at 
the binary collisions for proper numerical treatment of the classical equations of motion. In 
Sec. Ill scattering time signals are given as a function of the initial state of the hydrogen 
atom at the instant the positron starts coming to the hydrogen atom at a hxed distance. 
Based on observation of the dips in the scattering time signal, a code system is introduced 
for classifying orbits, and a surface of section (SOS) is also introduced for analyzing the 
phase space structure. In Sec IV we present our hndings based on the calculational result of 
the partition of SOS showing the phase space structure. Discussions are also given in that 
section. Finally, we conclude with a summary in Sec. V. 


II. EQUATIONS OF MOTION 


We consider the collinear p'^e~e~^ conhguration in which the proton and the positron are 
located on the opposite side of the electron. We further assume that the proton is inhnitely 
heavy and hxed at the origin. In atomic units, this system is described by the Hamiltonian 


H 


1 2 2 1 

-Pl - PiP2 +P2 - 

2 xi 


1 1 

- 1 - 

X2 Xi + X2 


( 1 ) 


where xi and X 2 are respectively the distances of the proton and the positron from the 
electron, and their conjuate momenta are denoted by pi and p 2 - 

For proper numerical treatment, the singularities at Xi = 0 or 0:2 = 0 are regularized by 
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introducing a fictitious time r and a transformed Hamiltonian K such as [21] 


dt = XiX2dT 

K = xiX2{H - E) 


( 2 ) 

(3) 


In addition, introducing regularized coordinates Qi, Q 2 , and momenta Pi, P 2 such that 


Xi = Ql, X2 = Ql, Pi = 


Pi 


P2 = 


2 Qi’ "" 2Q2 

one can obtain the equations of motion as follows 


(4) 


Ql — ~^PiQ2 ~ ~^QiQ2P2 
Q 2 = —^QiQ 2 Pi + -^P 2 Q\ 


Pi — -^Q2PiP2 + 2Qi — -P 2 Q 1 + “^QiQt 
P 2 = \PiQ 2 {QiP 2 - Pi) + 2 Q 2 + 2 QIQ 2 


E- 

E 


Q 


{Ql + Ql)^\ 
Ql 


iQl + Ql)\ 

where the dot notation is used for the hctitious-time derivative. 


(5) 

( 6 ) 

(7) 

( 8 ) 


III. SYMBOLIC CODE SYSTEM AND SURFACE OF SECTION (SOS) 

It has been found useful to examine scattering time signals as a starting step to systematic 
investigation of the classical dyamics of three-body problems [22l |23|. In Fig. a typical 
signal for scattering time is shown as a function of the angle variable (0) of the action- 
angle variable pair of the one-dimensional hydrogen atom at time f = 0. The energy of the 
hydrogen atom was chosen to be Pi = —1.2, and the total energy of the scattering system 
to be P = —1. The positron started inward at the distance of 100 from the proton, i.e. at 
xi + X 2 = 100 and we computed the trajectories until the positron reached xi+ X 2 = 1000 
after scattering. Note that the scattering time is not the hctitious time but the real time 
spent by the trajectory. 

The scattering time signasl show typical chaotic scattering signals with inhnitely many 
peaks and dips, and we focus on the dips as we did in two-electron atoms[23]. Note that 
two wings stretch from the primary dip located at 0 = 0^ ~ 12.61. as can be seen in Fig. 

The left wing stretches beyond 0°, and its end reaches a boundary (at 0 ~ 240°) of the 
region of chaotic signals; the right wing stretches to the other bounday at 0 ~ 150°. The 
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FIG. 1: (Color Online) The scattering time as a function of the angle variable ((/>) for the total 
energy E = —1 and the energy of the hydrogen atom Ei = —1.2. 

wings rise infinitely at the ends. And, by examining many other dips in the region of chaotic 
signals, we conhrmed that every dip has its own two wings, and every wing rises infinitely 
at its end. 

In Fig. 1^ two scattering trajectories are depicted with the initial conditions slight de¬ 
viating from the primary dip; (a) one with the initial condition (j) < 0^, and (b) the other 
with (f) > (pfi- From these plots, we can see that the trajectory with the initial condition 
corresponding to a dip would undergo the triple collision and furthermore the triple collision 
separates two qualitatively different scattering products, i.e. inelastic and exchange scatter¬ 
ing as in Figs. |^a) and (b) respectively. By examining trajectories with initial conditions 
around many other dips in the region of chaotic signals, we conhrmed that every dip is pro¬ 
duced by the trajectories ending in triple collision and every pair of wings stretching from a 
dip are branches to different scattering products. 
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FIG. 2: (Color Online) Scattering trajectories evolved with the initial conditions (a) (p = 12.58 < 4>d 
and (b) (p = 12.64 > pd- Solid (red) line for electron’s position xi as a function of time, and dotted 
(blue) line for positron’s position x\ + X 2 - Insets are magnifications around the instant when 
trajectories nearly approach the triple collision. 

Symbolic dynamics is a useful tool for a classification of orbits in nonintegrable systems 
[21112S]- From the observation of the above-mentioned separation of dynamical behavior 
by the triple collision, we naturally introduce a symbol code system by dehning a sym¬ 
bolic assignment such as Symbolic dynamics is a useful tool for a classihcation of orbits in 
nonintegrable systems [2H [25] . From the observation of the above-mentioned separation of 
dynamical behavior by the triple collision, we naturally introduce a symbol code system by 
dehning a symbolic assignment such as 
' 

1 if the trajectory undergoes a binary collision of the proton and electron 

2 if the trajectory undergoes a binary collision of the positron and electron 
c if the trajectory undergoes the triple collision 

For example, the scattering orbit in Fig. [^a) can be represented by a symbol sequece 
...111212111... while the orbit in Fig. [^b) by ...11121222... It is noticeable that the triple 
collision is a non-regularizable singularity and thus we cannot evolve trajectories any more 
after undergoing the triple collision mm. Thus we let triple collision orbits terminate with 
the symbol c. 

It is convenient to explore the phase space structure by using the Poincare section method 
[71I2SII2S]. It is noted that, by virtue of the energy scaling property of Coulomb systems, we 
only need to deal with the system for a hxed total energy E = —1 Based on the obvious 
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FIG. 3: (Color Online) Partition of SOS induced from classification of orbits by symbol sequences 
of length 2. The SOS consists of two sheets, xi = 0 (the right half-plane) and X 2 = 0 (the left 
half plane). Points represented by different symbol sequences are painted with different colors. 
Variables Q's and P’s are the regularized coordinates inroduced in Sec. II. 

fact that all orbits except the shortest CTCO undergo binary collisions {xi = 0 or 0:2 = 0) 
at least once, the set Xi = 0 or X 2 = 0 composed of two sheets of binary collisions is chosen 
as the SOS (see Fig. |^. The Poincare mapping is then naturally dehned as F{X) = V if a 
trajectory starts at X G SOS consecutively crosses at V G SOS. 

A partition of SOS can be induced by combining the Poincare mapping F and the symbol 
code system introduced above. For an orbit starting from X G SOS, a symbolic sequence of 
hnite or inhnite length can be given to it as 

S{X) = (ToCTi . . . CTfc . . . 
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where 


if F^{X) lies in the sheet Xi = 0 
if F^{X) lies in the sheet X 2 = 0 

if the orbit starting from F^~^{X) falls in triple collision 
without any more crossing of SOS. 

In Fig. 1^ we present a coarse-grained partition of SOS, where the whole SOS is parti¬ 
tioned to four regions of the sequences of length 2. The borders of the regions are the stable 
manifold of the triple collision which consists of orbits terminating with the symbol c. In 
the following section, we will present further hner partitions of SOS, which make it possible 
to hnd remarkable dynamical properties of the system. 

IV. RESULTS AND DISCUSSION 

Proceeding by three step from the Fig. we obtain a finer partition of SOS in which 
each region is classihed by a sequence of length 5 as shown in Fig. Comparing the two 
sheets (a) xi = 0 and (b) X 2 = 0, we can hnd that they have the same geometric structure. 
Actually we conhrmed this symmetry of geometric structure by examining the partition 
classihed by symbolic sequences of length 7, a plot of which is not given here for brevity. 
Thus we propose a conjecture that the geometric structure of the dynamics is invariant 
under exchange of the variables (xi,pi) with {x 2 ,P 2 )- This is a very surprising hnding since 
the masses of proton and positron are quite diherent and thus equations of motion (5)-(7) 
is not invariant under the exchange. 

Checking the symbol sequences representing the partitioned parts in Fig. we can see 
that there are no trajectories representing the symbol sequences 11122, 11221, 11222, 12211, 
11211, and their rehections 22211, 22112, 22111, 21122, 22122. By examining the hner 
partition up to symbol length 7, we conjecture that there is a hnite set of forbidden patterns 
of symbol sequences such as 1122, 2211, 11211, 22122. Note that there are no scattering 
orbits ...111222... or ...1112111...; the simplest scattering orbits are ...111212111... and 
...11121222... as can be seen in Fig. 

The partition of SOS is also useful for understanding the main features of the scattering 
time signals in Fig. Initial conditions evolve with time and pass through the SOS. A 
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FIG. 4: (Color Online) Partition of SOS by symbol sequences of length 5 repsented in two sheets 
(a) xi = 0 and (b) X 2 = 0. Colors are used for the same purpose as in Fig. 
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FIG. 5: (Color Online) Intersection of SOS and orbits with the initial conditions(/> G [5° : 15°] and 
4> G [150° : 240°]. The former is denoted by A; the latter by B. The segments A and B are plotted 
in the SOS partitioned by symbol sequences of (a) length 5 and (b) length 7. 

set of initial conditions at their first passage through the SOS forms an initial line segment 
in the SOS. And then the initial line segment proceeds in SOS by the Poincare map. In 
Fig.|ll we show images of two initial line segments obtained by applying the Poincare map 
several times. The segment A is the image of initial line segment formed from the interval 
0 G [5° : 15°] which includes the initial condition <f>d corresponding to the dip with the widest 
wings in Fig. the other (denoted by B) for the interval 0 G [150° : 240°] enclosing the 
initial conditions corresponding to the chaotic signals consisting of many peaks and dips in 
Fig. 0 It can be seen from Fig. |^a) that both the the line segments A and B are cut by 
the stable manifold of tripl collision. However, looking into the partition of SOS by symbol 
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sequences of length 7, we can find the number of intersection points of the segment B and 
the stable manifold of triple collision increases as the partition gets hner as in Fig. |^b). 
This is consistent with the feature of the scattering time signals such that there are inhnitely 
many dips in in the interval 0 G [150° : 240°]. 

From studies of three-body problems, it has been generally known that trajectories in 
captured state can exit to uncaptured state along only the orbits falling in triple collision; the 
closer the trajectory gets to triple collision, the larger the escaping velocity of the trajectory 
gets HISS]. In other words, the farther the trajectory gets from triple collision, the smaller 
the escaping velocity gets. This means that orbits with initial conditions far from dips, i.e. 
very close to the ends of wings escape the three-body interaction zone with zero energy of 
relative motion between the hydrogen atom (positronium) and electron (proton). Thus the 
ends of the wings rise inhnitely, producing peaks in scattering time signals as was mentioned 
in Sec. III. 

In addition to the role of triple collision as an exit from captured state, the other promi¬ 
nent role of triple collision is led by time-reversal symmetry of the present system as follow: 
scattering trajectories can enter into captured state along only the orbits starting from triple 
collision [221 [23]. Thus, in semiclassical physics, resonances in scattering cross sections can 
be described in terms of contributions from CTCOs, i.e. closed orbits starting and ending 
in the triple collision [3l jOj [20112^. It can be easily seen that all CTCOs lies in the collinear 
conhguration, and thus the present study of classical mechanics of collinear system 

is relevant to the huctuations in the scattering cross sections in the real 3-dimensional space. 
Therefore, as was done for two-electron atoms PEI, signihcant predictions can be made on 
the total scattering cross section of positron by the hydrogen atom with the total energies 
approaching the three-body breakup threshold from below as follows: (1) fluctuation in the 
cross section as a function of the total energy would be chaotic, (2) the amplitude of the 
fluctuation decays algebraically as the energy approaches the three-body breakup threshold, 
and (3) the Fourier transform of the fluctuation would show distinct peaks at the position 
of the actions of CTCOs. 

The statement (1) comes from the chaotic nature of the symbolic dynamics of the system: 
between any two ordinary admissible symbol sequences, we can hnd arbitrarily long symbol 
sequences with ending in triple collision. The statement (2) comes from the fact that the 
instabilities of relevant part of CTCOs approach to inhnity as the energy approaches the 
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FIG. 6: A phase portrait revealing an island of stability in the chaotic sea. The center of the tori 
is indicated by an arrow. 

breakup threshold, which is related to the irregularizability of the triple collision [SI [9]. The 
exponent for algebraic decay is determined by the stability of the relevant part of CTCOs, 
which will be obtained in future work. Finally, the statement (3) is based on closed orbit 
theory [2E], which is a general theory for atomic photoabsorption cross sections. However, 
the condition required in the theory is nothing but the restriction on the starting and ending 
point of the leading orbits to captured state: both the starting and ending point should be 
close to the center of atomic system. Thus it can be applied to scattering cross sections for 
the systems with CTCOs. The action values and relative instabilities of CTCOs are not 
calculated in the present paper, remaining as a future work. However, we can predicted that 
there are no peaks in Fourier transform of the fluctuation which correspond to the forbidden 
symbol sequences cll2c, cl22c, c221c, c211c, etc. This is quite important characteristics 
different from two-electron atoms, reflecting the phase space structure of the present system. 
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Looking into the partition of SOS, we get doubt if there are some islands around the 
point {Q 2 ,P 2 ) ~ (1.64,0.0) in Fig. In order to lift obscurity, we draw a phase portrait 
around that point. As can be seen in Fig. there is actually an island of stability consisting 
of tori in the chaotic sea. At the center of the island which is indicated by an arrow in Fig. 

there is a stable periodic orbit, which is found to be the shortest one with inhnite symbol 
sequence 12; all orbits in the island is represented by the same symbol sequence. This is 
another characteristics different from two-electron atoms of which collinear eZe subspace is 
known to be fully chaotic [3 [23. Since it has been accepted that there is no bound state 
of positron-hydrogen system in the real 3-dimensional world [TB] , the stability of the tori in 
the direction orthogonal to the collinear subspace should be investigated for semiclassical 
prediction of existence or absence of bound or long-lived states. And it will be also a good 
subject of future study to understand effects of the tori through dynamical tunneing in the 
scattering cross sections [29] . 

For complete study of positron-hydrogen system in future, we should mention on a scat¬ 
tering channel which is excluded in the present work, i.e. positron-electron annihilation. 
It is well known that the possible contribution of the positron-electron annihilation to the 
total scattering cross section is negligible except in the limit of zero positron energy HU. 
However,it may be necessary to consider positron-electron annihilation for total energies just 
below the three-body breakup threshold since the amplitude of fluctuating part in the total 
cross section would decay to zero as the total energy approaches the breakup threhold and 
thus would be affected even by small perturbations owing to the annihilation channel. 


V. CONCLUSION 

The classical dynamics of positron scattering by the hydrogen atom in the collinear con- 
hgurational subspace was studied for total energies below the three-body breakup threshold. 
A systematic analysis of the phase space structure of the system was given for the hrst time, 
resulting in important new hndings: (1) the system shows chaotic scattering, (2) the topo¬ 
logical structure of the phase space is invariant under exchange of the dynamical variables 
of proton with those of positron, (3) there is a hnite set of forbidden patterns of symbol 
sequences such as 1122, 2211, 11211, 22122, (4) there are tori forming a stable island around 
the shortest periodic orbit. And, based on the phase space structure and well-known general 
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dynamics near triple ollision.we explained main features of scattering time signals. 

Using the relevance of CTCOs to fluctuations in scattering cross sections owing to reso¬ 
nances, quantum manifestations of the classical dynamics were predicted on the total scat¬ 
tering cross section such that (1) fluctuation in the cross section as a function of the total 
energy would be chaotic, (2) the amplitude of the fluctuation decays algebraically as the 
energy approaches the three-body breakup threshold, and (3) the Fourier transform of the 
fluctuation would reflect the phase space structure through absence of peaks represented 
by the symbol sequences cll2c, cl22c, c221c, c211c, etc. including the forbidden symbol 
sequences. 

This work is considered as a hrst step to full understanding of classical, semiclassical and 
quantum dynamics of positron-hydrogen system, and we suggested several subjects to be 
done in future. In addition, it would be a kind of challenge to hnd a quantum mechanical 
menifestation of the structural invariance of the phase space under exchange of the dynamical 
variables of proton with those of positron. 

Acknowledgments 

This work was supported by the Kumoh National Institute of Technology under contract 
number 2011-104-015. 


[1] G. H. Wannier, Phys. Rev. 90, 817 (1953). 

[2] G. Tanner, K. Richter, and J. M. Rost, Rev. Mod. Phys. 72, 497 (2000). 

[3] M.-H. Lee, G. W. Byun, N. N. Choi, G. T. Tanner, Phys. Rev. A 81, 043419 (2010). 

[4] R. McGehee, Invent. Math. 27, 191 (1974). 

[5] J.-M. Rost, Phys. Rep. 297, 271 (1998). 

[6] H. Friedrich, Theoretical Atomic Physics (Springer, Heidelberg 2006), Chap. 4. 

[7] K. Richter, G. Tanner, and D. Wintgen, Phys. Rev. A 48, 4182 (1993). 

[8] Y. H. Jiang, R. Piittner, D. Delande, M. Martins, and G. Kaindl, Phys. Rev. A 78, 021401(R). 
(2008). 

[9] C. W. Byun, N. N. Choi, M.-H. Lee, and G. Tanner, Phys. Rev. Lett. 98, 113001 (2007). 


12 



[10] W. H. Cherry, Secondary electron emmision produced from surfaces by positron bomardment 
(Ph.D. Thesis, Princeton University 1958). 

[11] M. Charlton and J. W. Humberston, Positron Physics (Cambridge University Press, Cam¬ 
bridge 2001). 

[12] C. M. Surko, G. F. Gribakin, and S. J. Buckman, J. Phys. B 38, R57 (2005). 

[13] A .K. Bhatia, J. Atom. Mol. Cond. Nano Phys. 1, 45 (2014). 

[14] J.-M. Rost and E .J. Heller, Phys. Rev. A 49, R4289 (1994). 

[15] W. Ihra, J. H. Macek, F. Mota-Furtado, and P. F. O’Mahony, Phys. Rev. Lett. 78, 4027 
(1997). 

[16] A. S. Kadyrov, I. Bay, and A. T. Stelbovics, Phys. Rev. Lett. 98, 263202 (2007). 

[17] K. Jansen, S. J. Ward, J. Shertzer, J. H. Macek, Phys. Rev. A 79, 022704 (2009). 

[18] O. R. Spivack, J. Phys. IV France 9, Pr6-191 (1999). 

[19] Phys. Rev. A 77, 044502 (2008). 

[20] G. Tanner, N. N. Choi, M.-H. Lee, A. Czasch, and R. Doner, J. Phys. B 40, F157 (2007). 

[21] S. J. Arseth and K. Zare, Celest. Mech. 10, 185 (1974). 

[22] N. N. Choi, M.-H. Lee, and G. Tanner, Phys. Rev. Lett. 93, 054302 (2004). 

[23] M.-H. Lee, G. tanner, and N. N. Choi, Phys. Rev. E 71, 056208 (2005). 

[24] B.-L. Hao and W.-M. Zheng, Applied Symbolic Dynamics and Chaos (World Scientific, Singa¬ 
pore 1998). 

[25] Z.-Q. Bai, Y. Gu, and J.-M. Yuan, Physica D 118, 17 (1998). 

[26] M.-H. Lee, N. N. Choi, and G. Tanner, Phys. Rev. E 72, 066215 (2005). 

[27] M.-H. Lee and N. N. Ghoi, J. Kor. Phys. Soc. 56, 1799 (2010). 

[28] M. L. Du and J. B. Delos, Phys. Rev. A 38, 1896; ibid. 38, 1913 (1988). 

[29] S. Tomsovic and D. Ullmo, Phys. Rev. E 50, 145 (1994). 


13 



